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THE EFFECTIVE SPIN HAMILTONIAN AND PHASE SEPARATION 
INSTABILITY OF THE ALMOST HALF-FILLED HUBBARD MODEL AND 

NARROW-BAND s-f MODEL 
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Abstract 

The effective spin Hamiltonian is constructed in the framework of the almost half-filled Hubbard model on the Cayley tree by 
means of functional integral technique with the use of static approximation. The system in the ground state appears to be 
consisting of the ferromagnetic metallic domains and the antiferromagnetic insulating one sprovided that the concentration of 
excess electrons (or holes) does not exceed some critical value. The connection between the Hubbard model and the s-f model 
is stated. 



TO DESCRIBE the narrow-band magnetic semiconduct- 
drs, one uses in general two models: (1) the Hubbard 
model with the Hamiltonian 



gated to the spin density (see reviews 



ija 



(i) 



e\ (r)dr 



Q{^ (r)} 



where c. , Ci a are the creation and annihilation electron 

operators on the site i with the spin projection a, = 

J 



(3) 

where j3 1 = T is the temperature (K = 1, k = 1), 
Q{ei(r)} is the partition function of free electrons in 
the external magnetic field Eiir)- Then we use a com- 
mon "static approximation" |g, LZ[ , i.e. we replace £j (r) 



c\ a c ia , fj, is the chemical potential, u > \tij\, an average by its mean value over (0, (3) - £j. At t tj = 0, Q {£i (r)} 



number of electrons per site N e /N = 1 + c (c > 0) being 
close to unity; and (2) the s-f model 



can be calculated exactly. Considering formally the band 
energy term in equation ([1} as a perturbation we get 



T~t s -f — 
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+ H f 



(2) 



where Si is the /-spin operator, cr is the Pauli matri- 
ces, TLf is the Hamiltonian of direct /-spin interactions, 
/ is the exchange parameter, |/| ^> \Uj\, N e /N <C 1. 
The first model describes the situation when the conduc- 
tion electrons belong to the same band as the electrons 
forming the localized moments. On the other hand the 
Hamiltonian (U) is applied when these electrons are in 
the different bands (s and /). Nevertheless the models 
(fT|) and @ are equivalent under certain approximations 
as it will be proved further (Nagaev was the first who em- 
phasized this feature Therefore we restrict ourselves 
to consideration of the Hubbard model. 

When both c and \tij\/u are small the competition 



of tendencies to ferro- and antiferromagnetic orderings 



x [cosh (Pfi) + cosh (/?£,-)] exp (— /?<!> {ej}) , 
${£j} = -TTrln(J-tG), 



(4) 



where fi = [i — u/2, t, G are the operators with matrix 
elements: 



Gf (r,r') = (^n)e- <w "( T - T ') 



respectively, u n = irT (2n + 1), n = 0, ±1, ±2, 



> c 



(iu> n + J1) —A 



(5) 



takes place 
instability 




which leads to the phase separation 
_ . So far this effect has been treated qual- 
itatively (see also [H, Q). We shall construct here the 
effective spin Hamiltonian for the systems under consid- 
eration which will give the possibility to describe this 
instability as the first-order phase transition. 

We proceed with the representation of the Hubbard 
model partition function Z in terms of the functional 
integral over the fluctuating vector fields {ei (r)} conju- 



is the Green function of zeroth order approximation in 
tij . $ = $ {ej } is the functional of the spin interaction 
energy. When calculating the integral over the modulus 
of Sj, we can apply the saddle point method since u T, 



eorij (eo = u/3, \nj\ = 1). Then $ appears 



to be the functional only of unit vectors rij. The formal 
perturbation expansion of $ {n^ } with respect to Uj may 
be obtained in the usual way. It turns out that the true 
parameter of the expansion is \Uj \ / |erj — It is of order 



1 



of \tij \ /eo <C 1 at c = 0. But if c is small but finite, c » 
exp(— /3eo), we have ju = eo + Tlnc/(1 — c) in the zeroth 
order approximation with respect to , and therefore 
the perturbation expansion fails. In the case of c ^ we 
use the representation of functional ${n.,} through the 
exact Green function Qij (uj\X): 



1 f 

$ = / dujf {u)-jL)y j Im Aj , w + = uj + iO, 

Ai {uj) = Tr CT / dXy^tjjQjj (uj\X) = 
Jo 



y — Tr CT [(uj - eon, • a) Qu (uj[\) - 1] 



(0) 



in which /(a;) = (exp f3x+ 1) 1 and H A) satisfies the 
following equation 



(uj - e n; ■ <r) Q. tj {uj\\) - A/,, 0',- , (a;|A) 



(7) 



5iij there 



For the s-/ model with classical /-spins Sj 
is the correspondence with equations (6) and (7) if one 
replaces eo — > |J| 6 1 . This proves rigorously the equiva- 
lence of the narrow-band Hubbard model in the static 
approximation to the s-f model with classical /-spins 
and Hf = 0. 

Because the equation (|6]) contains the diagonal Green 
functions only we try to derive a closed equation for it. It 
is convenient to introduce the locator self-energy Ci (w|A) 
by the relation 



Qu 1 MA) - uj - 



eon; 



A MA) 



(8) 



The following expression for A (uj[\) was proposed in |8| 
(for the case of electrons moving in disordered media) 

3 

(At) 2 5> +W HA) (9) 



where 8 labels the nearest neighbours (the nearest neigh- 
bour approximation is assumed for iy). Equation ([9|) 
can be obtained from the second-order term of the per- 
turbation theory for A MA) by replacing Gj HA) — > 
Qjj {uj\X). It turns out to be exact on the Cayley tree 
[8J. Both in the ferromagnetic (FM) and the anti- 
ferromagnetic (AFM) cases • iij+s does not depend 
on i. Therefore we assume that the function cos6*j 
= (l/z)X)(5 n i ' n i+s {z is the coordination number) is 
a smooth function of site i. In this approximation and 
on account that \t\ /eo <C 1 the matrix equation ((§]) can 



be solved. The solution to equation (J9j) gives 

A, {uj) = 2 J — [{uj 2 - e 2 Q ) X, HA) - l] 
Xi (w|A) = 2 {uj 2 - e 2 + (uj 2 - e 2 ) 2 - 8z (Xt) 2 

xM + e^cos^)] 172 }" 1 (10) 

(the branch of square root is chosen by the condition 
Xi (w|A) — ► (uj 2 — el) 1 at t — > 0). Substituting equation 
(jTDJ) to equation taking into account the smallness of 
c, |t| /eo and integrating over A and uj we obtain 



e(x) 




157T 



const; 



1, a; > 0, 
0, x < 0, 



where 



= e — 2y/z \t\ cos ■ 



(1 — cos 6^ 



(11) 



(12) 



$ {n 3 } is desired free energy functional. It is necessary to 
write down the equation for the shifted chemical potential 
ju in addition to equations (JTTJ) and (TT2")) . It reads 



= (c{nj}) $ , c{nj} = 



3ttA^ 



E 



3/2 



|t| cos ; 



xGQi-Ei); 



i 

/n dn J ex p(-/ 3$ { n j}) 



(13) 



At T = the total energy of the system E = -T In Z + 
ptA e equals £ = min(${nj-} + Ncp,{rij}), ju{nj} being 
the chemical potential for the fixed spin configuration 
when c — > c {rij}. 

We minimize the energy putting cos = 1 in the do- 
main of relative volume x (FM phase) and cos If = in 
the domain of relative volume 1 — x (AFM phase). The 
optimal value of x is x = c/cq, 



Co 



±V /3 

2 V 3tt J e 



I|t| 



3/5 



(14) 



At c < c the two-phase state obtained has the en- 
ergy which is lower than that of any homogeneous state. 
Apart from numerical factors, equation (|14l) agrees with 
the result of the qualitative consideration [4] . It can be 
seen from equation (|13p that all excess electrons are in 
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FM domains. At c > cq the homogeneous FM ordering is 
favourable. Our Hamiltonian reduces to double-exchange 
Hamiltonian derived earlier in s-f exchange model (see, 
e.g. Q) provided that cq <C c <C 1. It is worthwhile 
to mention that the Cayley tree approximation gives the 
right value for the AFM indirect exchange parameter but 
overestimates the energy of FM ordering in \%fz times. 

The contribution of space derivatives of the function 
cos#i to the free energy functional has been calculated 
also. It appears that the thickness of the boundary be- 
tween FM and AFM phases is of order of lattice constant. 
This matter will be published elsewhere. 
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